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An automatic system to calculate two loop weak corrections to muon anomalous magnetic 
moment is discussed. Diagrams are classified into eight types, according to their topology. 
Adopting Civitanovic-Kinoshita representation of Feynman amplitude and using the 
i topological property, the renormalization is performed consistently by n-dimensional 

, regularization method. 
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1. Introduction 

, GRACE system 1 has been extended to produce one loop Feynman amplitude in au- 

tomatic way and the order a corrections to several 2 — > 2 reactions are calculated 2 . 

. The GRACE can also generate higher loop diagrams and correponding amplitudes. 

As a first application of the system to two-loop calculation, we consider the weak 

■ corrections to muon anomalous magnetic moment. 

D , QED corrections to muon (<?-2) are calculated up to eighth order and its valus is 3 

Aa i Q ) ED = (g-2)/2 = 1165846947(±46 ± 28) x 1(T 12 

The first error arises from theoretical uncertainty and the second one from a mea- 
surement uncertainty of the fine structure constant a. The next order corrections 
in QED are estimated as AaqJ D = (39 ± 10) x 10~ n . The vacuum polarization 
due to hadrons contributes Aahadron = (703 ± 19) x 10~ 10 . The weak interaction 
calculated at one loop level contribute Aa^ l e a ° k p = (195 ± 1) x 10 -11 . The two loop 
weak correction is estimated in the leading logarithmic approximation and it reaches 
about 22 % of the one loop result 4 . So, it is worth to get two-loop contributions 
from 1678 diagrams completely in automatic way. 



X 



2. Formalism 

We adopt the Civitanovic-Kinoshita 5 representation of the Feynman amplitude. 
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Combining the denominators by Feynman parameters, the general muon vertex at 
two loop level (in the case of 6 internal lines) is expressed as 

where, pj is a momentum flowing on the internal line j and expressed as 

pj = E v.uy.+Qj- ( 2 ) 



8=1,2 



The symbol f] s (j) represents (1,-1,0) according to the flow of loop momentum t s and 
qj is a external momentum flow on the line j. The F(D) represents a numerator 
and is expressed by several momenta and 7 matrices. By shifting the loop momenta 
to delete linear terms, the denominator becomes 

Denom. =^U st (£ s ■ £ t ) - V (3) 

s,t 

z i m ) ~ H z ibi ' + IrTn ^ ^^(U ■ 9j), (4) 



! ' " ' detU 

1,3 



where 



and 



U s , t = Y,ZjVs(j)VtU) (5) 



Bij = (0% U^st 1 ■ dctu = Bji- (6) 

These matrices depend only on the topology of the diagram. By changing the basis 
of the independent loop momenta, the matrix U is diagonalized and we rescale the 
loop momenta. We adopt n-dimensional regularization method to remove the Ultra 
Violet divergence. As the numerator can be reproduced by the differential-integral 
operator Dj, we first integrate over n-dimensional loop momenta. The operator 
(j represents the line number) is defined as follows. 



1 f°° 3 



The basic relation to pick up the numerator is 

1 

Hp*-™*) (pi-miy () 

The fundamental relations necessary for us are given by 
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DW j_ = W + (_^_\ (10) 



-'3V m ' V m ' V 2 detuy (to - 1) V*"- 1 



where 



The numerator can be reproduced by replacing p^ 1 with the operator Dj. The basic 
formula for our calculation becomes (except for coupling constant), 

r(6 - n) iff 



(detU)™/ 2 - ief- n 



r(5-n) Ff r(4-n) 



-2 (dctU)™/^ 1 (V - ief~ n 4 (detU)"/ 2 + 2 (V - ie)^ n 
Practical method to generate proper numerator is as follow. 
• replace p^ in numerator — > (cfc^ + Q' 



(12) 



c° terms -> Fq 



• c 2 terms -> Ff : replace fcf fcj -> (-1/(2 detU))^^ 

• c 4 terms — > : replace 

K^AK - (-V(2 dctU)) 2 +g li xg vl ,B' ik B' jt + g^g vX B' u B' jk ] 

In the above replacement, the contraction of Lorentz indices should be done in 
n-dimension. We introduce the following projection operator in 4-dimension to pick 
up the contribution to (g — 2).(<7:photon mom., (p — q/2), (p + q/2):muon mom.) 

Proj(M) = \(jt> + m)(m 7 „(p.p) - (m 2 + ^)p M )0 + \i + m) (13) 

By this projection, the F 2 term drops out. We replace F ,Fi with the projected 
one /0, /l. The factor (—1/2) is included in fl. 

3. Removal of the Ultra Violet Divergence 

The ultra violet (UV) divergence arises from the one loop sub diagram 6 . The 
over all divergence (F 2 term) is irrerevant in our case, because of the projection. 
It is necessary for order (a 2 ) charge renormalization. The diagrams are classified 
into 8 different types of topology. They are shown in the Fig.l. We will discuss two 
types of topology as examples. The vertex and self-energy type. (See Fig. 2.) 

3.1. Vertex type 

We divide the diagram into a sub diagram (S) and the remained (R). The 
correponding Feynman parameters are 

(z2,z 3 ,z e ) £ 5", (z l7 z 4 ,z 5 ) £ R (14) 
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The determinant (detU) has a semi-factorized form in the parameters. This implies 
that the singularity occurs when all the parameters in sub-diagram S tends to 0. 
In order to see this clearly, we reparametrize the Feynman parameters as follows. 

z 2 =x(l-y), z 3 =xy(l- z), z 6 = xyz (15) 

z\ = (1 — x)(l — u), Z4 = (1 — x)u(l — v), z 5 = (1 — x)uv (16) 

The measure of integration is changed into x 2 (l — x) 2 yu. When x approaches 
then the detU also tends to 0. This is the origin of divergence. It is important 
that the structure of determinant only depends on the topology of the diagram and 
the type of different topology is restricted in small number. So we can prepare 
the subtraction procedures beforehand, according to the type of topology. The UV 
divergence occurs from the second term. 

• the second term /l 

Setting e = 2 — n/2, the main part of the integrand becomes 

r^ ^^^^/V )^,,.,- (.7) 

where detU = x ■ u(x,y,z). (We must keep numerators up to 0(e).) We can see 
the existence of the singularity at e = when x — > 0, due to the factor x 6 ^ 1 . In 
order to extract the singularity, let us consider the following integral. 

/ = / x e - 1 {h (x)+£-h 1 {x))dx (18) 
Jo 

Partial integfration and the expansion in e give 

/=iM0) + M0)-£log*^M (19) 
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Using the above formula, the contribution from fl term turns out (as the integrand 
of dydzdudv{yu}) 

if'M = (i + 2^-2, E )fm + fm ( .^ gvm+ im 



e bV ' ' a J V(0) V(0) y b v " V(0) 



L 



l dx \u(x,y,z) 3 V(x) J 

in unit of (a/n) 2 /16. 

3.2. Self energy type 

Similarly, we divide the diagram into a sub diagram (S) and the remained (R) . 
The correponding Feynman parameters are 

(22,25) S 5", (21, 23, 24, 2 6 ) S i? (21) 

Reparametrization of Feynman parameters is done in the similar way. 

Z2=x(l—y), 25 = xy (22) 

21 = (1— x)(l— u), 23 = (1 — x)u{l — v), 24 = (1 — x)uv(l — w), z^ — (l~x)uvw 

(23) 

The measure is slightly different from the previous case, x(l — x) 3 u 2 v. The singu- 
larity occurs even in the first term /0. The treatment is the same as the fl in the 
previous case. The second term has the differnt structure. The main part of the 
integrand is 

To discuss the singularity, let us consider the following integral. 

1=1 x s - 2 {g (x)+e- gi {x))dx (25) 
Jo 

Partial integration and the expansion in e give 



/= 1%M + %^ + %M_ 5o(1) _ [\ og x^dx (26) 
e dx dx dx J dx 2 - 



By this formula, we can extract UV divergent part and the finite contribution. 

4. Counter terms 

There are 58 diagrms including one loop counter term. The contributions from 
these diagrams are easily obtained, by interpreting the counter term as the new 
coupling including (1/e). No divergence occurs in one loop integration, however, 
we must keep order e quantity in the integration to get final result. 
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vertex type self-energy type 

Figure 2: Assignment of Feynman Parameters 

5. System of Calculation 

We first prepare the topology dependent part. Corresponding to each topology, 
we can assign the flow of loop momenta and external momenta. Then we can 
calculate the important quatntity detU, Bij and the coefficients expressing the 
weight of the flow of momenta q and p in each internal line by using the conservation 
law of momentum at each vertex. These are the prepared files. The GRACE 
generates diagrams and correponding FORM source codes. By invoking FORM, 
the prepared files are included according to their topology and FORTRAN source 
codes giving the UV and finite part contibution to (g — 2) are generated. 
Numerical integration over Feynman parameter space is done by BASES 7 . For test 
run, we have checked several diagrams. We have reproduced the pure QED results 
by this system. The UV divergence is extracted systematically. The infrared (1R) 
singularity is regulated by introducing a tiny photon mass. The diagrams including 
self-energy type sub-diagram develop superficial IR singularity due to the fact that 
the x-derivative to get rid of UV part increases the power of denominator function. 
We need subtle treatment for IR singularity. 
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